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~qq ". Abstract 

Study of the fusion reactions between hydrogen isotopes in muonic molecules is the 

£NJ , first example of the accurate accounting of the nucleus charge screening by a muon in 

the fusion process. At LUNA installation the measurements of astrophysical reaction 

cross sections were extended down to collision energies of a few keV. Lhe screening 

by atomic electrons of the target became substantial. The possibility to look over 



< 

screening from unbound electrons is given by metal-hydrides used as targets in old 
reaction measurements. The classical Debye screening in plasma, applied to quasi- free 
electrons in metal, provides an explanation of unexpectedly large screening potentials 
found for some metals in the research through the Periodic table of elements. 
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1 Electron and muon screening of fusion reactions 

Fusion reaction cross sections decrease exponentially with decreasing of collision energy. It 
is convenient to isolate this dependency and characterize the reaction with the astrophysical 
function S(E) rather than with the reaction cross section o~(E): 

S(E)=Eo-(E)exp(27rr ] ) J (1) 
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where r\ is the Sommerfeld parameter, 2tci] = 31.29ZiZ 2 (m/ E) ' 5 , Z\ and Z 2 are charges 
of colliding nuclei, m is their reduced mass in amu, E is the cm. collision energy in 
keV (h = c = 1). Usually S(E) is a smooth function of energy. However, the success 
in the extrapolation of the measured S'-function to the energy of the astrophysical scale 
is not guaranteed: first, due to possible occurrence of bound or resonant states at low 
energies, second, due to the effect from atomic electrons of the beam or the target in collision 
experiments [lj. The latter effect leads to the increase of the cross section in comparison 
with the cross section of the interaction between "bare" nuclei which is expressed in terms 
of an enhancement factor 

f(E) = a cxp (£)/a barc (£), (2) 

where Obare and <T e xp are the cross sections for bare nuclei and the one measured in the 
presence of target and projectile atomic electrons. 

Electrons can be thought of as effectively lowering the Coulomb energy between the two 
colliding nuclei by a constant and energy-independent increment U e , the screening potential, 
which is usually treated as the acceleration of the incident particle by the electron cloud. 
In the adiabatic limit, the accelerating potential, £/ a d, is evaluated from the difference in 
atomic binding energies between the compound atom and the projectile plus the target 



atoms of the entrance channel. For reaction D(D,p)t the accelerating potential is equal 
to Ugd = 2 • 13.6 eV = 27.2 eV due to atomic electrons at the Bohr radius. The electron 
screening effect was observed for reactions d(d,p)t and d(d, n) 3 He with a molecular target 
in [2] . The extracted magnitude of the screening potential, 
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25 ± 5 eV, 



qualitatively agreed with the above estimate (Fig. 1). 

At lowest energies currently accessible in the laboratory, 1-10 keV, E 3> U e 
hancement factor behaves exponentially 

f{E) = a barc (£ + U e )/cr hsxe {E) « exp^r,U e /E) > 1 



(3) 



the en- 
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and even a small increase of the cross section at E/U e m 100 can lead to a wrong extrapo- 
lation if the screening correction is not properly introduced. 
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Figure 1: Comparison of the experimen- 
tal d(d,p)t S(E) -functions [2] with the 
bare-nuclei S , (£')-function from Ref. [3] 
(solid line) and estimate [4] (dashed line) 
in which electron screening is added to 
the bare-nuclei S- function using Eq. (4). 
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Figure 2: Measured astrophysical S(E)- 
function of d(d,p)t obtained with a 
deuterated Ta foil (T=-10°C) for atomic 
and molecular deuteron beams [15J. The 
dotted curve represents the bare S(E)- 
function, the solid one includes the effects 
of electron screening with U e =340 eV. 



In order to reduce the uncertainties behind extrapolation procedures in collision exper- 
iments, considerable efforts have been spent to push the experimental limits towards lower 
energies. LUNA collaboration measured the cross sections of the reactions 3 He(<i,p) 4 He 
and <i( 3 He,p) 4 He at the cm. energies, respectively, E = 5 to 60 keV and 10 to 40 keV 
and studied the screening effect [5]. The magnitude of screening potentials for both reac- 
tions, U e = 219 ± 7 and 109 ± 9 eV, appeared to be significantly larger than the values 
from atomic physics models, U e = 120 and 65 eV. This is not understood so far. The 
experimental analysis generally assumes that the bare-nuclei cross section, e.g., [3], derived 
by the extrapolation of data at higher energies, which are little affected by screening, is 
known. Authors [5 J conclude that, besides additional theoretical work on screening, a direct 
measurement of bare-nuclei 5'(£')-function would be desirable to clear up the situation. 



A possibility to verify the extrapolation is given by the muon catalyzed fusion. Muonic 
molecule is an ideal model system, which makes possible study of fusion reactions screened 
by a muon at practically zero energy [5J. The rates of nuclear fusion from the muonic 
molecule bound rotational-vibrational states (Jv), X"[ v , are expressed through the "bare" 
reaction constants Kl or the "bare" ^-functions 

\ J L V = K L ■ G{\ K L = S L (0)/(7t(am) 2L+1 ). (5) 

Here Sl is the S'-function, corresponding to the L-wave partial cross section, L is the orbital 
angular momentum of nuclei, Gj? is the density probability of nuclei penetration to the 
fusion region. This quantity, the "muon screening factor", was accurately calculated for 
muonuc molecules (see [7] and refs. therein). 

At PSI, radiative capture rates from the ground (00) state of pd\i muonic molecule were 
measured for doublet and quartet nuclear spin states [8]. At TUNL, doublet and quartet 
cross sections of reactions p(d, 7) 3 He and d(p, 7) 3 He with polarized nuclei were obtained for 
energies down to E cm = 23 keV and extrapolated to E = [9J. The values So(0) for the 
s-wave deduced from these experiments: So(0) = 0.109 ± 0.01 [9] and Sq(0) = 0.105 ± 0.01 
[5J (in units of eV b), are in an agreement, the extrapolation being thus verified. 

LUNA collaboration measured the d(p, 7) 3 He cross section from 22 down to 2.5 keV 
cm. energy, well below the solar Gamow peak [ID], where the cross section falls to values 
about picobarns. Their result for the total ((s + p)-wave) extrapolated S'-function 5(0) = 
So + Si = 0.216 ±0.010 is somewhat larger than that deduced from [9J S(0) = 0.166 ±0.014. 
The difference is within the uncertainty due to electron screening effect: the estimated 
enhancement of the S-factor is about 6% at 2.5 keV (cm.) and it increases to 20% for 
interaction energies around 1 keV, but screening potential could not be extracted from the 
data HDJ. 

Characteristic features of charge-symmetric reactions d(d, n) 3 He and d(d,p) 3 H: ratio 
of channels in p- and s-waves and p-wave contribution to the total cross section in both 
channels were studied via muon catalysis. Ratio R(n/p) of neutron and proton fusion 
yields in the dd\i molecule was measured as a function of deuterium mixture temperature 
[HJ. In conditions of resonant ddfi formation (p-wave fusion) the value R p (n/p) = 1.42 ±0.03 
was obtained. This ratio R(n/p) decreases to R s = 1.06 ± 0.04, when non-resonant ddfi 
formation and s-wave fusion become prevailing. The p-wave reaction constants for mirror 
channels were extracted from the fusion rates determined in [Sj. All results agree with the 
R- matrix analysis of (id-fusion in-flight data [3] , proving the efficiency of /iCF approach. 

Properties of dd reaction were studied in [IT] at a beam energy as low as 5 keV and 
target deuterons embedded in the metallic environment. The angular distributions of both 
2 H(<i, p)t and 2 H(d, n) 3 He reactions were measured. They are expected to be anisotropic 
even at the lowest energies due to a relatively large p-wave component in the entrance 
channel. The angular distribution is symmetric with respect to 90° and can be parametrized 
as A{9) = 1 ± a2cos 2 (0) (neglecting an L=2 contribution). The energy dependence of 
the parameter 02 for both reactions recommended in [TJ] was confirmed at extremely low 
energies. This agrees with the parametrization of [TJ], with the experimental results [13] 
and [2], and with the fiCF result [8]. 

The ftCF method might be helpful for understanding electron screening in 3 He(<i,p) 4 He 
and <i( 3 He,p) 4 He nuclear reactions discussed above. The fusion rate in muonic molecule 
<i/i 3 He was calculated using advanced methods of solving three-body Colomb problem 



Its measurement would provide the "bare" reaction constant, which is necessary for ade- 
quate fitting of the in-flight data [5] . 

2 Electron screening effects in solids 

The application of metal hydrides as targets provides a unique possibility to investigate 
the influence of free (unbound) electrons on fusion cross sections in beam-target collisions. 
During past decade several groups have performed a series of experiments bombarding 
deuterated metals with low energy deuterons. The electron screening effect in d(d,p)t 
reaction for the target produced via implantation of low-energy deuterons unexpectedly 
appeared to be about one order of magnitude larger than the value (3) from a gas target 
experiment [2]. As an example, Fig. 2 from [T5] shows the astrophysical S- function of 
the d(d,p)t reaction measured in deuterated Ta. The measured points clearly show the 
exponential increase at low energies. The derived screening potential for Ta is U e = 340 ± 
24 eV, which is significantly larger than the theoretical predictions based on the adiabatic 
approach. 

Early studies of deutron-induced in-flight reactions on deuterated targets were initiated 
because of the interest to cold fusion experiments. An anomalous enhancement of the fusion 
yield was reported [16] for deuterated Pd (U e = 250 ±15 eV) and a deuterated Au/Pd/PdO 
multilayer (U e = 601 ± 23 eV), while deuterated Ti and Au exhibited a normal (gaseous) 
enhancement: U e — 36 ± 11 and 23±11 eV, respectively, u 

In [TTJ thick target yields of the fusion reactions 2 H(d,p) 3 H. and 2 H(<i, n) 3 He were mea- 
sured on deuterons implanted in three different metal targets (Al, Zr and Ta) for beam 
energies ranging from 5 to 60 keV (Fig. 3) . The values of the screening energy also demon- 
strated a clear target material dependence U sc = (190 ± 15) eV for Al, (297±8) eV for Zr 
and (322±15) eV for Ta. @ 

Study by LUNA collaboration [15] of the reaction d(d,p)t for deuterated Ta tested 
surprising results on large screening energy of [TT] with somewhat different experimental 
approaches and confirmed them. This stimulated systematic research through the Periodic 
table and results for several metals, insulators, and semiconductors were reported in [17] 
and completed in [18], [19] (Fig. 4). 

A comparison of the U e values with the Periodic table indicates a common attribute: 
for each group of the Periodic table, the corresponding U e values are either low (gaseous) 
as for groups 3 and 4 and the lanthanides, or high as for groups 2, 5 to 12, and 15 (Fig. 5). 
Group 14 is an apparent exception to this characteristic: the metals Sn and Pb have a high 
U e value, while the semiconductors C, Si, and Ge have a low U e value indicating that high 
U e values are a feature of metals. A similar situation is found for group 13: B = insulator, 
Al and Tl = metals. The indication is supported further by the insulators BeO, Al 2 03, and 



-"■The screening effects found for deuterated Pd-containing targets could not justify cold fusion results. 

2 At the same time, no target material dependence of the angular distribution was observed in |llj . 
Thus the screening effect noticeable for beam energies below 25 keV should be equal for both entrance 
wave function contributions L = and L = 1 and should not change the relative intensities of both 
contributions compared to the gas target experiments. The ratio of the total yields of mirror reactions 
appeared close to unit for all three metal targets at all energies down to 5 keV, in the agreement with 
low-energy extrapolations and ^CF result. Cold fusion observations are evidently contradicting to this 
firmly established feature. 
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Figure 3: Enhancement of the thick- 
target yields for three different deuter- 
ated metals Al, Zr and Ta from [TTJ (2£ = 
deuteron lab. energy). Extracted values 
of the screening potential are shown. 
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Figure 4: S^-E 1 ) -function of d(d,p)t reac- 
tion on deuterated Cu, Nd, Hf, Pt (E in 
c.m.s.) [19]. The dotted curve shows the 
bare S(E) -function, the solid one repro- 
duces the exponential enhancement with 
electron screening potential U e . 



Ca02, as well as by deuterated metals M having an observed small stoichiometric x value 
(M X D) of the order of one or smaller and thus representing also insulators (e.g. group 4 of 
the Periodic table and the lanthanides) . In summary, a large screening effect is observed 
in all metals except in the noble metals Cu, Ag, and Au. 

Various aspects of the metals were discussed to explain the data. In particular, effects 
from deuteron thermal and vibrational motion in the lattice was found to be marginal 
[20J. Other mechanisms considered in [T5], [T7] - channeling, diffusion, conductivity, crystal 
structure, electron configuration, and Fermi shuttle acceleration - were unsuccessful either. 

The observation that large enhancements have been found in deuterated metals but not 
in insulators has suggested a possible explanation based on effects of the plasma of electrons 
in the metal [TS|, [IS]. This model with classical quasi-free valence electrons predicts an 
electron screening distance around positive singly charged ions (deuterons in the lattice) of 
the order of the Debye length 
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(A;T/(4Tre 2 n e //Pa)) 1/2 , 



(6) 



where T is the temperature of free electrons, n e ff is the effective number of valence electrons 
per metallic atom, p a is the atomic density. At room temperature and at n e ff w 1 one 
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Figure 5: Periodic table showing the studied elements, where those with low U e values 
(U e < 100 eV, small effect) are lightly shadowed and those with high U e values (U e > 100 
eV, large effect) are heavily shadowed [19] . 



obtains a radius Rd, which is about a factor of 10 smaller than the Bohr radius of a 
hydrogen atom. With the Coulomb energy between two deuterons at Rd set equal to U e , 
one obtains U e = 300 eV, reproducing the correct size of the screening potential U e . The 
acceleration mechanism of the incident ions leading to the high observed U e values is, within 
this simple model, the Debye electron cloud at rather small radius Rd- The n e ff values were 
compared with those deduced from the known Hall coefficient: ChoII = (en e //(ifa//)p a ) _1 
|21j : within 2 standard deviations the two quantities agreed for all metals. 

Some consequences of the Debye screening model were tested experimentally: 

• Debye screening energy Ud is proportional to target nucleus charge. Tested for pure 
metals Z = 3 (Li), 4 (Be), 23 (V), 71 (Lu) [22], [23]. 

• The Debye screening energy Ud is proportional to the nuclear charge of the incident 
ion. This prediction was verified [19] in the d( 3 He,p) A He studies in metals (Zi = 
2): with a typical value of U e = 300 eV for the dd fusion in metals at T = 290 K 
one would expect for d( 3 He,p) 4 He the Debye value Ud = ZiU e (dd) = 2 x 300 eV = 
600 eV, U e = [Tad + U D = 120 + 600 eV = 720 eV, f/ ad being the adiabatic screening 
potential; observation [19] gives U e = 680 ± 60 eV. 

• The electron screening is the effect of the entrance channel. The two mirror reactions 
measured with neutrons and protons in the exit channel [UJ, [23] demonstrate that 
the electron screening is not influenced by the the charged ejectiles of the exit channel. 

A critical test of the classical Debye model would be the predicted screening potential 
dependence on the temperature: U e oc T -1 / 2 . The electron screening in the d(d,p)t reaction 
has been studied for the deuterated metal Pt at a sample temperature T = 20°C -r340°C 
and for Co at T = 20°C and 200°C 0. The astrophysical S(£)-function obtained at 
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T = 20°C and 300°C for Pt is shown in Fig. 6. The U e (T) values for Pt are plotted in 
Fig. 7 together with the expected dependence U e oc T -1 / 2 (dotted curve). All data show a 
decrease of the screening, i.e. the U e value, with increasing temperature. 
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Figure 6: The S (E) -function of d(d,p)t 
for Pt at T = 20°C and 300°C, with 
the deduced solubilities y. The curves 
through the data points include the bare 
S(E) factor and the electron screening 
with the given U e values |2i 
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Figure 7: The observed values U e (T) 
for Pt are shown as a function of sam- 
ple temperature T. The dotted curve 
represents the Debye model prediction, 
the solid curve includes the observed T- 
dependence of the Hall coefficient |2l 



It was noted in [19] that the metals with a high hydrogen solubility, of the order of one, 
have small screening potentials. Contrary, for the metals with high U e values, the solubilities 
are quite small. Actual solubilities were not always available at room temperature, but those 
observed in [19] were consistent with known results [25] . 

The temperature dependence of this correlation was further studied. The electron 
screening effect for deuterated metals with high solubility, such as Ti (group 4), was mea- 
sured at temperatures from 10 to 200°C. The deduced solubility y(T) shows a sizable de- 
crease with increasing temperature. Above 50°C the solubility reaches values below 1 and 
an enhanced screening becomes observable. Finally, all metals of groups 3 and 4 and the 
lanthanides have been studied at T=200°C Similarly as for Ti, all these metals exhibited 
a large reduction in solubility and thus showed a large screening. U Screening potentials U e 
and solubilities y for 49 metals in total were measured in [T9], [21] (see also |26|). 



3 In the insulator C the solubility decreased to 0.15 (T=200°C) from 0.35 (T = 20°C), but no enhanced 
screening was observed, as expected for an insulator with n e ff = 0. 



3 Discussion 

All data on the enhanced electron screening in deuterated metals can be explained quantita- 
tively by the Debye model applied to the quasi-free metallic electrons. So far there is no idea 
why the simple Debye model appears to work so well in this case. It was argued [19] that 
most of the conduction electrons are frozen by quantum effects and only electrons close to 
the Fermi energy {Ep) actually should contribute to screening, with n e ff(T) ~ kT/E F oc T 
and thus there should be essentially no temperature dependence for U e . However, this 
argument applies only to insulators and semiconductors with a finite energy gap, while for 
metals there is no energy gap and the Fermi energy lies within the conduction band. 

Possible many-body effects characteristic for dense matter resulting from internuclear 
correlation processes can distort the picture; but, due to trapping deuterons by the crystal 
lattice, these effects can be neglected. One should refer to the series of papers by Ichimaru 
et al. on nuclear fusion in various types of dense plasmas (see the review [27]). In particular, 
different effects in metal deuterides were considered. They found that the electrons both 
in metallic d-band and hydrogen-induced s-band contribute to screening and it is more 
efficient in Pd than in Ti. However, the magnitude of the screening potential later observed 
in [25] remains unexplained, as well as its temperature dependence as in |24j . 

Without theoretical grounds, the Debye model can be treated as a powerful parametriza- 
tion of the data with an excellent predictive power. 

Authors [26] noted the correlation between the screening energy and the deuteron den- 
sity in the host during the bombardment. They relate the density to the diffusivity, or 
mobility of D + ions in the host; large mobility results in small density. The inverse of 
density is related to a concept they call the "fluidity" of the target deuterons. They pro- 
pose that high fluidity in the host might be responsible for the enhanced screening. The 
fluid deuterons and conduction electrons might behave like a plasma in the host and both 
electrons and positive ions reduce the Coulomb repulsion. A model of |28j is cited in this 
context. It is argued that the screening energy due to the fluid deuterons can be one order 
of magnitude larger than that due to the electrons because of the difference between Boson 
(deuteron) and Fermion (electron) statistics. Thus, another dynamic screening mechanism 
during the deuteron bombardment and penetration into the host is suggested, wherein the 
fluidity of deuterons must play a decisive role. 

To explain the observed fusion rate enhancement an approach based on the quantum- 
uncertainty dispersive effect between energy and momentum proposed by Galitskij and 
Yakimets [29] has also been applied. Authors [30] showed that the quantum-tail effect 
produces an important increment of the rate enhancement at very low energies, however, 
this mechanism cannot reproduce the observed experimental rate for deuterium. 

We conclude with the remark that the influence of the environment in various nuclear 
processes is one of the interesting subjects because of its interdisciplinary nature involving 
nuclear physics, condensed matter physics, and material science. Thus, low-energy nuclear 
reactions at far below the Coulomb barrier should be explored more in various conditions, 
experimentally as well as theoretically. Finding the method to describe the screening may 
eventually help understanding processes in stellar or terrestrial plasma. 
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